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O ■ Abstract 

^: 

pL^ . We analyze the influence of decaying sterile neutrinos with the masses in the 

^r) ' range 1 - 140 MeV on the primordial Helium-4 abundance, explicitly solving 

the Boltzmann equations for all particle species, taking into account neutrino 
flavour oscillations, and paying special attention to systematic uncertainties. 
' We show that the Helium abundance depends only on the sterile neutrino 

lifetime and not on the way the active-sterile mixing is distributed between 
^■i flavours, and derive an upper bound on the lifetime. We also demonstrate that 

the recent results of Izotov & Thuan [TJ, who find 2a higher than predicted 
by the standard primordial nucleosynthesis value of Helium-4 abundance, are 
consistent with the presence in the plasma of sterile neutrinos with the life- 
time 0.01 - 2 seconds. The decay of these particles perturbs the spectra of 
(decoupled) neutrinos and heats photons, changing the ratio of neutrino to 
qq ■ photon energy density, that can be interpreted as extra neutrino species at the 

recombination epoch. 

<N 

o: 

1 Introduction 

The characteristic feature of the physical processes in the early Universe is a pe- 
culiar interplay of gravity and microscopic physics. Gravity introduces the Hubble 
time parameter th that indicates the timescale on which the global properties of the 
Universe (geometry, temperature, etc.) change significantly. The Hubble time is de- 
termined solely by the energy density of the matter filling the space. The microscopic 
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matter constituents, particles, are involved in the interaction processes, that are be- 
lieved to be described fundamentally by three known forces — electromagnetic, weak 
and strong. As long as the timescale r of any given microscopic physical process is 
much smaller than r#, the expansion can be neglected on that timescale. If time r is 
enough to establish thermal equilibrium between the particles, then the equilibrium 
it maintained in the course of the Universe expansion, while r <C r# holds. When 
this inequality ceases to hold, the state of equilibrium is lost. The main reason for 
that is that interparticle distances become larger, while the corresponding densities 
become lower, hence interactions are less likely to occur. 

In this paper we are consider the formation of light nuclei in the primordial envi- 
ronment - Big Bang nucleosynthesis (BBN). All three fundamental interactions are 
important in the course of this phenomenon, all playing different roles. Charged par- 
ticles together with photons are subject to electromagnetic forces and the equilibra- 
tion timescale of corresponding processes is tiny with respect to the expansion time. 
Therefore the particles are kept in thermal equilibrium at the common temperature 
T. Due to expansion the temperature is decreasing with time. The equilibration time 
of the weak interactions changes abruptly so that at T > fewMeV weakly-interacting 
neutral particles (neutrinos and neutrons) stay in equilibrium, while at lower tem- 
peratures they fall out of it (freeze out). Weak interactions keep the difference of 
neutrons and protons, n n — n p , conserved, as in a sample reaction n + v e — > p + e + . 
At high temperatures the rates of this and similar processes are high enough to 
maintain chemical equilibrium, that is the neutron-to-proton conversion exhibits the 
same finite intensity as the opposite process. Chemical and thermal equilibria are 
interconnected for the weak interactions, so that they are lost simultaneously, at 
the moment indicated above. Finally, the strong interactions are responsible for the 
production of nuclei comprising more than one nucleon. The most important fusion 
reaction for the formation of the first nucleus, deuteron, n + p — )■ D, releases energy 
of at least the binding energy of deuteron Ed ~ 2.2 MeV, and proceeds effectively in 
dense primordial medium. As long as the temperature is higher than Ed, however, 
energetic photons collide with deuteron and lead to destruction of the latter. There- 
fore, the production of D is suppressed. As baryon density is much lower than the 
density of photons [2j, there are many photons with energies much higher than the 
temperature that collide with deuterons and hence postpone the production of the 
net deuteron density until the temperature T ~ 80 keV, that is much lower than the 
binding energy. Created D serves as a fuel for the formation of 3 He, 4 He and other 
nuclides. 

Although the time of the element production and the moment of the departure 
from the chemical p — n equilibrium are well-separated in time, the former process 
is very sensitive to the latter. Firstly, the details of the departure set the ratio 
of the neutron to proton densities and secondly, the time elapsed between the two 
moments determines the fraction of neutrons that have decayed since then (recalling 
that neutron is an unstable particle). 
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The seminal ideas of the primordial synthesis of light elements were first outlined 
in the so-called a/?7 paper, j_3j, published in the late 1940s. Since then the theory 
of Big Bang nucleosynthesis has evolved and its main predictions were confirmed, 
making it a well-developed model from both theoretical and observational points of 
view. A lot of reviews of the standard BBN scenario and its implication for particle 
physics models exist (see e.g. [1JEJE])- 

The predictions of the primordial nucleosynthesis can change once one replaces 
the Standard Model of particle physics underlying the processes considered so far by 
some of its "beyond the Standard Model" (BSM) extensions. Therefore the BBN 
plays the role of a benchmark for testing physical models. 

In this paper we investigate the influence of sterile neutrinos on primordial nu- 
cleosynthesis. Sterile neutrinos are hypothetical massive super-weakly-interacting 
particles (see e.g. [7J |8] for reviews), as opposed to their weakly-interacting coun- 
terparts - ordinary Standard Model neutrinos v e , z/ M , u T , that are called " 'active" in 
this context. Sterile neutrinos carry no charges with respect to the Standard Model 
gauge groups (hence the name), but via their quadratic mixing to active neutrinos 
they effectively participate in weak reactions and at energies much below the mass 
of the jy-boson their interaction can be described by the analog of the Fermi theory 
with the Fermi coupling constant Gf replaced by Gf x where the active- sterile 
mixing angle i?„ C 1 (see Fig. [p. Here a is a flavour index, a = e,/i, r, indicating 
that sterile neutrino can mix differently with neutrinos of different flavours. Massive 
sterile neutrinos can decay, but due to their feeble interaction strength their life- 
time can be of order seconds (even for masses as large as MeV). The decay products 
of the sterile neutrinos are injected into the primordial environment, increasing its 
temperature and shifting the chemical equilibrium. 

In this work we concentrate on sterile neutrinos with the masses in the MeV range, 
motivated by the recent observations |H] QUI [HI H2J E33] that particles with such masses 
can be responsible simultaneously for neutrino oscillations and generation of baryon 
and lepton asymmetry of the Universe and can influence the subsequent generation 
of dark matter [14J. The corresponding model has been dubbed z/MSM (Neutrino 
Minimal Standard Model, see [7J for review). 

Several works had previously considered the influence of MeV-scale particles on 
primordial nucleosynthesis. Compared to the Refs. [151 US] this paper accounts for 
the neutrino flavour oscillations in the plasma and employs more accurate strategy of 
solving Boltzmann equations, which results in the revision of the bounds of [T5"| [16] 
(see Section H] for detailed comparison). The authors of [17] developed a new code 
that can perform treatment of active and sterile neutrinos with arbitrary distribution 
functions, non-zero lepton asymmetry, etc. However, as of time of writing this code 
has not been made publicly available and the Ref. [17] did not derive bounds on 
sterile neutrino parameters. The work [18] concentrated on the bounds that cosmic 
microwave background measurements could provide on decaying sterile neutrinos 
with the masses 100 — 500 MeV, leaving BBN analysis for the future work. A number 
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(a) Quadratic mixing N v e 
of sterile neutrino to v e 
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(b) Decay of sterile neutrino N — > 
v e v a v a through neutral current interac- 
tions 



(c) Fermi-like interaction with the "ef- 
fective" Fermi constant i? e x Gp for the 
process in the panel (b). 



Figure 1: Fermi-like super- weak interactions of sterile neutrino 

of other works ([121 EDI 1211 1221 123]) analyzed the influence of decaying MeV particles 
on BBN. We compare with them in the corresponding parts of the paper. 

The paper is organized as follows. We explain the modifications of the standard 
BBN computations due to the presence of sterile neutrinos in the plasma and describe 
our numerical procedure in Sec.EJ Our results are summarized in Sec. [31 We conclude 
in Sec. HI Appendixes KHCl provide the details of our numerical procedure. 



2 Big Bang Nucleosynthesis with sterile neutrinos 

The section below summarizes our setup for the BBN analysis with decaying particles. 
Our notations and conventions closely follow the series of works [161 EU 122] ■ 

We will be interested only in the tree-level Fermi interactions of sterile neutrinos 
with the primordial plasma. In this case the interaction is fully determined by the 
squares of their mixing angles. We will consider one Majorana particle with 4 degrees 
of freedonQand three active- sterile mixing angles l d 2 a . Matrix elements of interactions 
of sterile neutrinos with the Standard Model particles are summarized in Appendix IBI 
(Tables [3] - 14 on page 23p . 



1 This number corresponds to g s — 2 of additional chiral singlets (i.e. "neutrino-like" species). 
Actual number of degrees of freedom is of course twice larger: dof = 2 x g s = 4, because every 
chiral fermion has 2 different helicity states. 
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We consider in this work only sterile neutrinos with the masses in the range 
1 MeV < M s < rj 140 MeV. For heavier particles, two-particle decay channels 
appear (e.g. iV — > 7i$v a) 7r ± e =F ) and our procedure of solving Boltzmann equations 
(described below) should be significantly modified. The lower bound is determined 
by the following considerations. The sterile neutrino lifetime t s is 



r- 1 = T S = [(1 + ~g\ + gl){$l + <%) + (! + £ + g 2 M] 



96vr 3 
6.9 sec -1 



10 MeV 



1.6^ + 1.13(^2 + ^2 



where 0w is the Weinberg's angle and gn = sin 2 ^ fh 0.23 , <7z = | + siia 2 9w, 
<3l = — |+sin 2 6^11 From this expression one sees sterile neutrinos lighter than 2 MeV 
have lifetime of at least several hundred seconds even for very large mixing angles 
•Q ~ 1 - much longer than the Hubble time at temperatures ~ 1 MeV. For smaller 
mixing angles lifetimes are even larger. Such longevous particles would increase the 
energy density of the Universe at the time of the freeze-out of weak reactions and 
therefore would effectively contribute to the number of relativistic neutrino species, 
N cS . The current bound puts N eS = 3.74+^ ± 0.06(syst) at 2a [lj (where the 
systematic error is due to the different values of neutron lifetime between the average 
value from Particle Data group, [27] and the recent measurement of |2B]). Therefore, 
one can derive robust bounds on the properties of sterile neutrinos with mass below 
MeV without solving a system of Boltzmann equations. Additionally, in the z/MSM 
the successful baryogenesis is possible only for the masses of sterile neutrinos above 
few MeV [HI [T3] . Therefore we restrict our analysis to the region of masses higher 
than 1 MeV. 



2.1 Expanding Universe and distributions of particles 

We consider expansion of the homogeneous and isotropic Universe with the flat 
Friedmann-Robertson- Walker metric in the form ds 2 = dt 2 — a^dx 2 , where a = a(t) 
is a time-dependent scale factor, whose evolution is described by the Friedmann 
equation 

a / 8ttG n 

H =a = y— p > (2) 

with the quantity on the left-hand side being the Hubble expansion rate, reciprocal 
to the expansion timescale th discussed above. The total energy density p is the sum 
of all the energy densities present in the medium, and Gn is the Newton's constant. 
The energy density together with the total pressure density p satisfy the "energy 

2 The expression (JTJ) is for Majorana particle. For Dirac particle the lifetime would be twice 
longer. 
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conservation" law 

a^ + 3(p + p) = 0. (3) 
da 

At the temperatures of interest the dominant components of the plasma are photons 
7, electrons and positrons e ± ) three flavours of active neutrinos (v e , v T ) and ster- 
ile neutrinosJl Working with the particle kinematics in the expanding Universe it is 
convenient to use conformal momentum y instead of the usual physical momentum p. 
The two are related through y = pa. The quantitative description of the plasma pop- 
ulation is provided by the distribution functions f a , that are the numbers of particles 
a per "unit cell" of the phase space d 3 p d 3 x = (27r) 3 . At keV-MeV temperatures the 
medium is homogeneous and the distribution functions are independent of spatial 
coordinates of particles, and due to isotropy they do not depend on the direction 
of the particle momentum. That simplifies the description of their evolution and 
therefore 

df fdf TT df\ df(t,y) 

i s \i - H "i) = ^ (4) 

holds. The goal is to find the distribution functions of all relevant particles and 
to use them to compute the energy density and pressure as a function of time and 
scale-factor, closing the system of Eqs. ((21)-© via 



fi(y,a)-J—dy (5) 
Ei{y) 



Here the summation goes over all plasma particles, mi is the number of degrees 
of freedom and mass of z-th particle respectively, Ei(y) = Ei(y)a = \Jy 2 + mfa 2 . 

If interaction rate of the particles is much faster than the Hubble expansion rate, 
their distribution functions are given by either the Bose-Einstein, or the Fermi-Dirac 
distributions. This is the case for photons, electrons and positrons — that are kept 
in equilibrium due to intensive electromagnetic interactions 

1 1 

A = e E/T _ l ' f e = e E/T + I ■ ( 6 ) 

Their contribution to the energy and pressure in Eqs. ©, © is hence determined 
by the single parameter - temperature. However, to describe the contributions of the 
other particles one has to solve kinetic equations involving them (see Sees. I2.3H2.4I 
below) . 



2.2 Baryonic matter 

The contribution of the baryonic matter to the evolution of the hot plasma of rela- 
tivists species is proportional to the so-called baryon-to-photon ratio r] B = n B /n^. 

3 Muons may appear in plasma from the decays of the sterile neutrinos with M s > 106 MeV. See 
Sec. HiU for details. 
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The measurements of relic radiation [2] yield T]b = (6.19 ± 0.15) x 10~ 10 . One can 
see that baryons are present in negligible amount, and do not influence the dynamics 
of the remaining medium. This allows to analyze our problem in two steps. At 
step I we omit baryonic species and study how the temperature of the plasma, the 
expansion factor and neutrino distributions evolve in time from temperatures of the 
order of 100 MeV, when sterile neutrinos typically start to go out of equilibrium^ 
down to Tpin ~ 10 keV when nuclear fusion reactions have ended. At step II we use 
these results to determine the outcome of the nuclear reaction network against the 
background of evolving electromagnetic plasma (Sec. 12. 5p . 



2.3 Active neutrinos at MeV temperatures 

Weak interactions are not able to maintain the thermal equilibrium of active neu- 
trinos with the plasma during all the expansion period we consider. A simple com- 
parison of the weak collision rate G 2 F T 5 and H(T) tells that neutrino maintain their 
equilibrium with the rest of the plasma down to temperatures T dec ~ few MeV. The 
process of neutrinos going out of equilibrium is usually referred to as neutrino decou- 
pling. Throughout the paper we assume that no large lepton asymmetry is present 
so that the number of neutrinos is equal to the number of antineutrinosH At tem- 
peratures higher than T^ ec the distribution is therefore given by the Fermi-Dirac one, 
while at lower temperatures we have to solve the set of three Boltzmann equations 

a = e, /j,t (7) 



dt 



1 a ■ 



The details of the interactions, such as particle collisions, are encoded in the so-called 
collision terms I a . The terms are explicitly [33] 



/~ = ^E i S\M\ 2 F[f}{2irY5\p m -p out )X[> ^ 



a ~ 2E a ^ J UJV ; KFm Fout; -Ll (2tt) 3 2^ v ; 

in, out i=2 

The sum runs over all the possible initial states "in" involving v a (represented by a 
particle set u a , 2, 3, . . . , K) and the final states "out" (K + l, . . . , Q). Matrix ele- 
ment Ai corresponds to the probability of the transition "in" -"out" to occur and the 
delta-function ensures the conservation of 4- momentum p; n = p ut- Symmetrization 
factor S is equal to 1, except of the transitions involving identical particles either in 
initial or in a final state. Relevant matrix elements together with the symmetriza- 
tion factors are listed in Appendix [Bj The interaction rates are dependent on the 
population of the medium, and the functional F[f] describes this. In case when all 



4 The order of the temperature depends on the mixing angle and can be even larger. 

5 For the previous studies of the BBN outcomes with the large lepton asymmetry present see 



e.g. [29, 3H33 [El 132] . 
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the incoming and outgoing particles are fermions, 



F[f] = (1 " UJ ••• (1 - f K )f K+ i ...f,y A,... f K {l - f K +i) ••• (1 - f Q ). (9) 



When some of particles are bosons, one has to replace (1 — /#) by (1 + /#) for every 
bosonic particle i?. 

A simple estimate (see Appendix [Cj demonstrates that the rates of transitions 
between neutrinos of different flavours are much faster than weak reactions. We argue 
that this phenomenon can be approximately described by the following modification 
of the Boltzmann equations 



Summation is carried out over three active flavours and expressions for Pp a are listed 
in Appendix 0(Eqs. [23]). 

2.4 The impact of sterile neutrinos 

As already mentioned, sterile neutrinos interact much more feebly than active neutri- 
nos do. Nevertheless, at some high temperature sterile neutrinos may enter thermal 
equilibrium. Whether this happens or not depends on the thermal history of the 
Universe before the onset of the synthesis^ Even if they were in thermal equilibrium 
at early times, sterile neutrinos then necessarily decouple at temperatures higher 
than those of active neutrino decoupling. If sterile neutrinos were light and stable 
(or very long-lived), they would be relativistic and propagate freely in the medium, 
yielding N e g « 3 + together with active neutrinos. However, sterile neutrinos 
decay into active neutrinos and other particles. The energies of the decay products 
may be very different from the typical energies of plasma particles. For particles 
that equilibrate quickly, (such as electrons or photons) this "injection" results in the 
fast redistribution of the energy between all particles in equilibrium and effectively 
the process looks like a temperature increase (more precisely, it just slows down the 
cooling of the Universe). But for particles that either are not in equilibrium or are 
about to fall out of it, such as active neutrinos at few MeV, the "injection" modi- 
fies the form of their spectra. The other mass-induced effect is that sterile neutrinos 
may switch from the relativistic regime (when their average momentum is larger than 
mass), that is established at large temperatures, to the non-relativistic one, due to 
the gravitational redshift. 

For the quantitative description of sterile neutrino dynamics we utilize the Boltz- 
mann equation similar to ([7]), replacing active neutrino everywhere therein by sterile 

6 For example in the ^MSM model at early times (T ^ 100 GeV) initial densities of sterile 
neutrinos are negligible |34) . Then the neutrinos come into equilibrium at temperature T + (typically 
T+ = 10 100 GeV) and freeze-out at temperatures T_ - 0.5 - 5 GeV [12]. 




(10) 
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neutrino 

Reactions contributing to the right-hand side together with their probabilities are 
listed in Tables |3]- 14 on page 23| of Appendix [B] Note that we neglect the processes 
involving baryonic particles. However, they become important for temperatures near 
the QCD crossover temperature Tqqd — 200MeV, when their density is not negligible 
anymore. More scattering channels of sterile neutrino would appear and their proper 
account is involved. However it seems to be reasonable to assert that the only 
modification the account will bring is to lower the decoupling temperature of sterile 
neutrinos. 

Oscillation phenomenon does not affect significantly sterile neutrinos and there- 
fore Boltzmann equation in its original form f ill I) is still valid, contrary to what 
we have found out for active neutrinos. An argument in favor of this statement is 
explained in Appendix ICl 

When sterile neutrino is heavier than muon, the former particle can appear in the 
decay iV — > n~ + e + + v e . However, the branching fraction of this decay mode does 
not even reach a percent for masses of sterile neutrino we consider (see e.g. [26J ) . 
Therefore we can neglect influence of both muons and other particles, appearing in 
the decay. 

As a result we have six equations (J2]), ([3]), ffTUj) . and ffTTj) describing primordial 
plasma at temperatures of interest. These equations contain six unknowns - scale 
factor a(t), temperature T{t) and four neutrino distribution functions/^ and f^. 
The system of equations is therefore closed and we have solved it numerically at the 

STEP I. 



2.5 Course of nuclear reactions 

Outcome of the nuclear reaction chain may be found numerically as well. For the 
Standard BBN model one of the earlier attempts was made with the code written by 
L. Kawano [35J. However, this program in its original form is inappropriate for the 
account of the BSM physics we consider in this paper, so we had to modify it. The 
resulting code became an essential tool of step II in our approach. Presence of sterile 
neutrinos alters the standard dynamics of the temperature and the expansion rate as 
well as the rates of weak interactions involving neutrons and protons. Knowing this 
from the step I, we used them as an input. The computations of nuclide evolution 
started from temperatures of several MeV, when the chemical equilibrium ceases to 
hold, up to temperatures T Fin . 
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2.6 Adopted values of Helium abundance 



The observables of the BBN are concentrations, or abundances, of light nuclides 
dispersed in the cosmos. The most relevant abundance in our problem is that of 
4 He as it is sensitive to the expansion rate of the Universe at MeV temperatures and 
neutrino distribution functions. The presence of sterile neutrinos in plasma typically 
increases the concentration of 4 He. On the contrary, deuteron is sensitive mostly 
to the value of baryon-to-photon ratio which is not significantly affected by sterile 
neutrino decays and interactions. In this work we adjust the initial value of baryon- 
to-photon ratio 7] at the beginning of the computation so that by T Fin ~ 10 keV 
it is equal to the value given by cosmic microwave background measurements [2J. 
Therefore the abundance of deuterons is essentially fixed in our framework and lies 
within the observational errors, not carrying the information about sterile neutrinos. 

The actual values of 4 He abundance are expressed traditionally in terms of mass 
fraction Y p , that is the fraction of total cosmological baryon mass stored in form of 
this element. Standard Model predicts the value of this abundance 

Y p SBBN = 0.2454 ± 0.0026 (syst.) (r n = 885.7 sec) 

y^sBBN = 2437 ± Q.0026 (syst.) (r„ = 877.7 sec) ' ' 

(see Appendix IA. 1 1 and Tabled] for details about systematic uncertainties in comput- 
ing this value. Notice that difference in Y p due to different values of neutron lifetime 
r n is smaller than our adopted systematic uncertainty. On r n see below). 

There are two main methods of experimental determination of primordial Helium 
abundance. The first one is related to the studies of low-metallicity astrophysical 
environments and extrapolating them to zero metallicity case. The Y p measurements 
are know to be dominated by systematic uncertainties. Therefore we adopt the Y p 
values from the two most recent studies, Refs. [H [36] that have slightly different 
implications. For recent discussion of various systematic uncertainties in 4 He deter- 
mination, see [3Tj . 

In Ref. P the value Y p = 0.2565 ± 0.0010(stat.) ± 0.0050(syst.) was obtained. 
Therefore, the 2a intervals that we adopt in our studies ar^| 

Y p = 0.2495 - 0.2635 (Ref. 0Q, 2a interval) (13) 

One notices that this result is more than 2a away from the Standard Model BBN 
predicted value of Y p , Eq. (fl2|) . 

Based on the different studies, the authors of [36] obtain^ Yp = 0.2574±0.0036(stat.)± 
0.0040(syst.). As a result 2a intervals become: 

Y p = 0.2462 - 0.2686 (Ref. [36], 2a interval) (14) 
7 We add the systematic errors linearly 

8 We use the average value over metallicities, (Y p ) (Eq. (8.2) of [35]) and add the difference 
between (Y p ) and its extrapolation to zero metallicity (Eq. (8.1)) as a systematic error. 
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Figure 2: Bounds (at 2a level) on sterile neutrino lifetime as a function of their 
mass for various measurements of Y p (summarized in Section [2TB]) . All results are for 
mixing of sterile neutrino with electron flavour only. 

(this values of Y p coincides with the Standard BBN one, f ll2p . at about la level). 

Second method of determination of Helium abundance is based on the CMB 
measurements. This method determines truly pristine value of Y p , not prone to the 
systematics of astrophysical methods. However currently its uncertainties are still 
much larger than of the other method. The present measurements put it at 

Y p = 0.225 - 0.401 (Refs. [381 12], 2a interval) (15) 

again consistent with the Standard Model BBN at 1.5a 

Finally, we mention another important uncertainty originating from the particle 
physics. There are basically two different measurements of neutron lifetime r n that 
are at tension with each other. Particle Data Group [27] provides r n = 885.7±0.8 sec, 
while measurements performed by Serebrov et al. [28] result in r n = 878.5 ± 0.8 sec. 
We employ both results and explore the differences they lead to in what follows. 

3 Results 

In this Section we present our results: the bounds on sterile neutrino lifetime as a 
function of their masses and mixing angles. As discussed in the previous Section, 
there are several systematic uncertainties in the determination of the 4 He abundance 
and therefore the results will depend on the adopted values of Y p (together with the 
neutron lifetime, r n ). We summarize these systematic effects below. 

We start with comparing the upper bounds on sterile neutrino lifetime for different 
values of Y p adopted (see Section 123]) . The left panel of Fig. [2] shows that the bounds 
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Figure 3: Upper bound for sterile neutrino lifetime for different mixing patterns: 
mixing with z/ e -only (red dashed line), z^-only (green dashed-dotted line) and equal 
mixing with v e and flavours (black solid line). Both bounds are derived by as- 
suming lifetime of the neutron r n from |28j . 

from the two recent works [1, 36J are quite similar (the difference is of order 30%) 
while the CMB bound is a factor ~ 5 weaker. 

The lower bound on Y p from the recent work of pQ is above the Standard BBN 
value ( Il2p at ~ 2a level (see however [37]). The presence of sterile neutrinos in 
plasma of course relaxes this tension and therefore at 2a the adopted values of Y p 
(Eq. [T3|) provide both upper and lower bounds on sterile neutrino lifetime. This is 
shown in Fig. |2j right panel. At 3a level the measurements of |T] are consistent with 
Standard BBN and the lower bound disappears. 

Another systematic uncertainty (the lifetime of neutron, r n ) is negligible. Indeed, 
the relative difference between sterile neutrino lifetimes were found to be of the 
order of 5% for two choices of r n - from [28] and from [27] (taking the same Y p 
bound from [36]). 

Next we investigate the dependence of the resulting bounds on the mixing pat- 
terns of sterile neutrinos. Naively, one would expect that sterile neutrinos mixing 
"only with v e " and "only with 2/ T " should produce different modifications of Y p . How- 
ever, due to the fast neutrino oscillations at the BBN epoch (see Appendix [Cj this 
is not the case and different mixing patterns give essentially the same results (see 

Fig. ED. 

4 Discussion 

In this work we considered the influence of decaying particles with the masses 1-140 
MeV on the primordial abundance of 4 He. Such particles appear in many cosmologi- 
cal scenarios HH EHl [231 [221 EU ISl SOI H21 131 SI]- Particularly, we concentrated 
on the properties of sterile neutrinos and derived constraints on their lifetime imposed 
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Figure 4: Comparison with the previous results of 



by the present measurements of primordial Helium abundance Y p . Sterile neutrinos 
are super-weakly interacting particles, quadratically mixed with the active flavors. 

We analyzed the case of one Majorana sterile neutrino with 4 degrees of freedom 
(if sterile neutrinos were kept in thermal equilibrium it would be equivalent to g s = 2 
species of active neutrinos). Due to fast neutrino oscillations at the BBN epoch, the 
bounds that we obtained are essentially independent on the mixing patterns as Fig. [3] 
demonstrates. 

In the paper [16] a similar model was considered with one Dirac sterile neutrino. 
Dirac sterile neutrino has the same 4 degrees of freedom and influences primordial 
plasma in the same way (if it has the same spectrum, lifetime and mixing pattern). 
However, in [T6] effect of active-neutrino oscillations was not accounted for. To 
provide corresponding analysis we wrote code that solves more accurate Boltzmann 
equations describing kinetics of neutrino than what were used in [16] . We compare 
our results with the previous bounds [151, HB] in Fig. SJ Our results are broadly 
consistent with the previous ones, however the differences can be as large as a factor 
of 2.5 for some masses. 

It is interesting to consider the parameters of sterile neutrinos, compatible with 
the lower bound, derived from the Helium abundance of [1] and compute the ratio of 
the total energy, deposited into the active neutrino species to that of photons (the 
quantity, proportional to effective number of neutrino degrees of freedom 



N, 



cff 



120 p Ve + pu^ + Pv T 



7tt 2 



rp4 



(16) 



at temperatures ~ 10 keV (when all sterile neutrinos had decayed). The results are 
plotted in Fig. [5j We see that N eS between 2.5 and 7 are possible for different mixing 
angles and masses, which can explain a larger than 3 values of iV e fj, reported recently 
in several CMB observations (see e.g. [3E1 S21 S3], but also [44]). Decaying sterile 
neutrinos with the masses 100 — 500 MeV and lifetimes from seconds to minutes and 
their influence of N e s and entropy production have been recently considered in [18] 
(see also [45]) where it was demonstrated that they can lead to N cS ^ 3 and can 
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Figure 5: Effective number of neutrino species (the ratio of neutrino to photon 
temperatures at T ~ fewkeV) as a result of decay of sterile neutrino with parameters, 
corresponding to the lower bound in Fig. [2j 

therefore be probed with the CMB measurements. Our results demonstrate that 
in the region 100 — 140 MeV where we overlap with the parameter space, studied 
in [18], the primordial nucleosynthesis restricts the lifetime of sterile neutrinos to be 
well below 1 sec (see Fig. |2j left panel). 

Finally, it is interesting to compare the upper bound on sterile neutrino lifetime, 
derived in this paper with the lower bounds that come from direct experimental 
searches for sterile neutrinos (see [26], HHJ HI])- We leave this for the future work. 
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Code 


Y p for r n from PDG [27] 


Y p for r n from [28J 


(Modified) Kawano code [35J 


0.2454 


0.2437 


PArthENoPE code [49] 


0.2480 


0.2462 


Difference 


-0.0026 


- 0.0025 



Table 1: Values of Helium abundance Y p in the Standard Model BBN (SBBN) and 
their uncertainties. 



A Tests of the numerical approach 

The Section below summarizes our comparison with previous work that analyze the 
influence of the MeV particles on primordial nucleosynthesis. Throughout this Sec- 
tion, we normalize scale factor by imposing condition aT = 1 at the initial moment. 
Conformal momentum is y = pa with the same normalization of the scale factor. 

A.l Standard Model BBN 

First we considered the nucleosynthesis in Universe filled with the Standard Model 
particles only. We compute the actual non- equilibrium form of the active neutrino 
spectra during their decoupling. Our results are compared with those of [25], 124] . 
Fig. |6] shows the evolution of the quantity aT as a function of temperature. It is 
identical to the Fig. 1 in Ref . [25] . Figure [7J shows relative distortion of neutrino 
spectra Sf Ue / f eq , Sf v / f eq , where fl q = (e v + is the equilibrium Fermi-Dirac dis- 
tribution and 5f u — f v — fl q describes the deviation from the equilibrium shape. One 
can see good agreement between the results and that the deviations from equilibrium 
are subtle. Also we compared asymptotic values of ratio aT at low temperatures. In 
[24] value 1.3991 is listed, while we derived 1.3996. 

The resulting Y p is summarized in Table [T] for different values of neutron lifetime 
r n . We also provide a comparison of the modified version of the Kawano code [35J that 
we adopted for computing nuclear reactions with a newer code, PArthENoPE |49j. 
According to [I] the systematic uncertainty of Kawano code is approximately +0.0017, 
so that actual value of Y p is larger by this value than what the code predicts. By 
comparing the results of PArthENoPE and the modified KAWANO code we find it 
to be 0.0026. 

A. 2 Test of energy conservation 

If all weak reactions involving electrons and positrons are turned off, neutrinos de- 
couple from the rest of plasma. Then the energy conservation law ([3]) holds separately 

fl http ://sec.bitp.kiev.ua| 
lc http : / / virgo . org . ua 
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Scale factor a [MeV '] Scale factor a [MeV '] 

Figure 7: Left: Relative distortion of v e spectra Sf Ue / f eq for conformal momenta 
y = 3, 5, 7 (from bottom to up). In each pair of curves the lower one corresponds to 
our result and the upper one corresponds to that of figure 3a from [23]. Right: The 
same, but for muon neutrino 
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for the neutrino component and for the remaining particles. In approximation of zero 
mass of electron we obtain 

similar to Eq. ( TT9|) . As a corollary, product aT is conserved. On the other hand, 
our code solves the equation ([3]) involving all medium components simultaneously. 
And it turns out that the relation ( ITT]) is not a trivial consequence of the numerical 
computation. Therefore the check of the conservation serves as a test of the code. 
We considered separately scattering and decay processes involving neutrinos and 
observed conservation of aT with precision of order 0.2%. 



A. 3 Heavy sterile Dirac neutrino 

Next we have tested model with one sterile Dirac neutrino iV with mass M s = 
33.9 MeV, mixed with v T [15J. This neutrino was assumed to be in thermal equi- 
librium with plasma at T > 50 MeV. To simplify the problem, the authors of [15] 
used the Boltzmann equilibrium statistics for active species in collision integral for a 
sterile neutrino. 

Being in equilibrium the sterile neutrino spectrum becomes more and more non- 
relativistic with time due to the redshift. Therefore the ratio p s /M s n s of the energy 
density p s to the mass times number density n s should approach 1 at lower tempera- 
tures. We have recomputed the evolution of the system using our code, without the 
Boltzmann approximation. Fig. [8] shows the comparison of the results with those of 
[TB] for sterile neutrino lifetime t s = 0.3 sec. Both results coincide till T « 5 MeV 
and after that moment ratio p s /M s n s of [T5] stops decreasing, while our numerical 
result shows the expected behaviour — the ratio continues to decrease, approaching 
1. 



A. 4 Massive v T 

Next we considered a model with the massive tau neutrino [191 EI]- Fig- El presents 
relative deviation of the energy densities of massless neutrinos 5p v / p CCL produced by 
our code and plotted in [T9"] . p eq = 7 ' K 1 ^ > is the equilibrium energy density of one 
neutrino specie, and 5p u = p u — pl q . In Fig. [9] distortion of electron neutrino spectrum 
y 2 6~f Ue / feq is depicted. Here one observes good agreement between the results. 



A. 5 Late reheating model 

To test the treatment of MeV decaying particles, we considered the low-reheating 
models with the reheating temperature of several MeV [221 E3]- In [22] heavy non- 
relativistic particles were considered, that dominated the energy density of the Uni- 
verse once and then decayed into electrons, positrons or photons (so that decay 
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Figure 8: Ratio p s /n s M s as a function of scale factor for M s = 33.9 MeV sterile 
neutrino. The upper curve is the result of Ref. [15], the lower curve is the present 
work. 
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Figure 9: Left: Relative deviation from its equilibrium value of v e energy density 
dpvjpeq in a model where tau neutrino is massive. Our numerical result is compared 
with [19J. Right: Spectrum distortion y 2 8f Ue /f e<l for the same model. In both panels 
M Vt = 0, 3, 7, 20 MeV from bottom to top. 
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Figure 10: Effective number of neutrino species N e g depending on decay width of 
heavy non-relativistic particles. Comparison of the results of this work and Refs. [22] 
and [23] . 

products are quickly thermalized) . The most important output is the effective num- 
ber of active neutrino species N e g (defined in Eq. f|T6|) ). Dependence of this quantity 
on decay width of heavy particle is presented in Fig. [TDJ We have noticed some 
difference between the results of cited papers and those of our code. We believe 
that this is due to the different approximations made. For example, in both works 
[221123] the scattering processes involving only neutrinos were not taken into account, 
approximation of Boltzmann statistics was used throughout and electron mass was 
neglected. We checked that the account of finite electron mass gives a gain of 5% to 
the N e s for r = 0.1s, while the account of scatterings involving only neutrinos gives 
rise of 1%. 

A. 6 Instant thermalization of decay products 

Next we considered a model with two heavy Majorana sterile neutrinos, similar to the 
z/MSM. However, we assumed that for any mass of sterile neutrino it can decay only 
via channels listed in Table H] of Appendix [Bl It is not a natural assumption, because 
usually sterile neutrinos heavier than pion decay dominantly into states containing 
mesons |26j. Also we approximated sterile neutrino spectrum as a non-relativistic 
one, while all the other particles are relativistic and in equilibrium all the time. In 
this case the system may be adequately described by the kinetic equation 




s 



(18) 
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Scale factor a [MeV -1 ] Scale factor a [MeV -1 ] 

Figure 11: Left: Evolution of aT for the model of Sec. IA.6I Right: p s / Psm- I n 
both panels semi-analytical calculations are compared with those of our code for 
three parameter sets: sterile neutrino mass M s = 580 MeV with lifetime r = lsec; 
M s = 1030 MeV with r = O.lsec; M s = lOOMeV, r = 0.5sec. 



together with the Friedmann equations (J2H3)). The latter of these equations can be 
rewritten as 

d(aT) _ 30aT s p s 
dt ~ 43tt 2 T 3 1 j 

T s is the decay width of sterile neutrino, p s is the energy density of sterile neutrinos, 
and we have used expression for the energy and pressure densities of relativistic 
species p re i = 3p re i = 437r 2 T 4 /120. 

In Figs. dUthe evolution of quantities aT and p s /p rc i is compared between the 
results of our code and the semi-analytic integration of Eqs. (!T8"]) - (fl9]) for three 
different sets of masses and lifetimes. One can see very good agreement between 
these results, maximum relative deviation is 1%. 



B Tree-level matrix elements 

In this Appendix we summarize the matrix elements we used for computing the 
collision integrals in Boltzmann equation. The squares of the matrix elements for 
Standard Model particles only are listed in Table El while the squares of the matrix 
elements of processes with sterile neutrinos are summarized in Tables [31 HI In these 
expressions, averaging over helicities of incoming particles and summation over those 
of outgoing products is assumed. The reactions are considered for two cases. In the 
first one sterile neutrino is a right-chiral Majorana neutrino that has 2 helicity degrees 
of freedom. That is actually the case in our problem, where we have two neutrinos 
f this kind. The other case corresponds to sterile neutrino of Dirac nature. Dirac 
fermions have both right- and left-chiral components, hence yielding 4 degrees of 
freedom in total. Expressions listed in Tables El 0] are applicable for both cases of 
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Process (1 + 2-^3 + 4) 


5 


SG F 2 \M\ 2 


77 -t- Ua — ^ 77 -1- IJq 

Vol \ u p 7 v a \ u p 


1 


32(pi -p 2 )(P3 -P4) 






v n A- Vr — \ Z7^ A- Vr 

u OL 1 " p ' ** OL 1 ** p. 


1 


32(pi -p 4 )(p 2 - Pa) 






1J A- 7/ — ^ 7/ -1- 7/ 


1/2 


64(pi -p 2 )(P3 -P4) 






Vs. A- v~ — \ V~ A~ Vs. 

OL ' OL ' OL I OL 


1 


128 (pi -p 4 )(p 2 -P3 


) 




7/ -4- 77 — V 77 a A- IJn 

Vol \ v a 7 i/ft ~r i/fj 


1 


32(pi -p 4 )(p 2 "Ps) 






V* A- Z7„ — V A~ f~ 


1 


128[s 2 (p!-p 4 )(p 2 - 


Ps) + 








0j?(Pi • Pa)(P2 • p 4 ) + 9l9r™1{P\ 


P2JJ 


v„ 4- c~ — y v„ 4- p~ 


1 


128[g 2 L ( Pl -p 2 )(p 3 - 


Pa) + 








9r(Pi -Pa)(P2 -Ps) 


-9L9Rml(pi 


P3JJ 


77„ -1- P - * - — ^ Z/„ -4- P~^~ 


1 


128[g|(pi.p 4 )(p2- 


Ps) + 








9 R (pi ■ P2)(P3 ■ Pa) 


-9L9Rml(pi 


P3JJ 


v, i \ A- Z7, i \ — >> A- f~ 

Vfl{T) \ V^ T j 7 O C- 


1 


128[g 2 ( Pl -p 4 )(p 2 - 


Ps)+ 








9r(Pi -Ps)(P2 -P4) 


+ 9L9Rm 2 e(Pi 


•P 2 )] 


^(r) + e~ ->■ i/ M ( T ) + e" 


1 


128[^( Pl -p 2 )(p3- 


Pa) + 








0h(Pi -P4)(P2 - Pa) 


- 9L9Rm 2 e(Pi 


■J*)] 




1 


128[^(p 1 -p 4 )(p 2 - 


Ps) + 








9r(Pi ■ P2)(P3 ■ Pa) 


- 9L9Rm 2 e(Pi 


•P3)] 



Table 2: Squared matrix elements for weak processes involving active species only. 
5 is the symmetrization factor; a,(3 = e,fi,r. In all processes we take a ^ (3. The 
results coincide with those of Ref . [25J . 

the neutrino nature. Moreover, to complete the list of possible tree-level reactions, 
one has to consider charge-conjugated channels and take into account that Dirac 
particle is distinct from its antiparticle, while Majorana neutrino is not. 

Throughout this Section we use the notations gn = sin 2 ^ , g L = 1/2 + sin 2 0\y> 
gL = —1/2 + sin^w, where 6 W is the Weinberg angle so that sin 2 6 W m 0.23. Our 
results coincide with [251 E] • 

C Neutrino oscillations 

The active neutrinos of different flavours v^v^ v T are related to the mass eigen-state 
basis U\, V2, V3 via a non-diagonal Pontecorvo-Maki-Nakagava-Sakata (PMNS) matrix 
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Process (1 + 2-43 + 4) 


5 


SG F 2 \M\ 2 


u s + v$ -4 V a + Up 


1 


6ZV a {Pi p 2 )\P3 Pa) 




Uc 4- ur — > i/„ 4- z/« 


1 


32^(pi-p 4 )(p2-p 3 ) 




I/ s + I/ a -4 u a + u a 


1/2 


0^iV a {Pl p 2 )\J>3 Pa) 




Uc 4- z/„ — > z/^ + u„ 


1 


128^( Pl .p 4 )(p 2 .p 3 ) 




U s + U a ^f Up + Up 


1 


6ZU a\Pl PA)\P2 P3) 




u s + P e — >■ e + + e~ 


1 


128# 2 e [g 2 (p 1 -p 4 )(p 2 -p 3 )+ 
g R (Pi ■ P3)(P2 ■ Pa) + 9L9Rml{p\ 


■P2)} 


+ e~ — )■ u e + e~ 


1 


128^ 2 e [g 2 L (p 1 -p 2 )(p 3 -p 4 )+ 
9r(Pi ■ Pa)(P2 ■ P3) ~ 9L9Rm 2 e {j)\ 


■Ps)] 


z/ s + e + — )• z/ e + e + 


1 


128$ 2 e [g 2 ( Pl -p 4 )(p 2 -p 3 )+ 
9r(Pi ■ P2){P3 ■ Pa) - 9L9Rm 2 e (pi 


■Ps)} 


v s + u^ t) -4 e+ + e~ 


1 


^K(rM(Pl -Pa){P2-P3) + 

9r(Pi ■ Ps)(P2 ■ Pa) + 9L9Rm 2 e(Pi 


■P2)] 


z/ s + e~ -4 z/ m(t) + e~ 


1 


^K(rM(Pl ■ P2)(P3 ■ P4) + 

9r(Pi ■ Pa)(P2 ■ Pz) ~ gL9Rm 2 e (Pi 


■Ps)] 


z/ s + e + -4 i/ m(t ) + e + 


1 


128$l {T) [g 2 ( Pl ■ p 4 )(p 2 ■ p 3 )+ 
g R (pi ■ P2)(ps ■ Pa) - hgRm 2 e (pi 


■Ps)] 



Table 3: Squared matrix elements for scatterings of sterile neutrinos N. Here S is 
the symmetrization factor; a, (3 — e, /j,,t; a ^ (3. d a is the squared mixing angle 
of sterile neutrino with u a . The results are applicable for one right-chiral Majorana 
neutrino as well as for one Dirac neutrino, for details see text. 
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Process (1 -» 2 + 3 + 4) 


5 


SG F 2 \M\ 2 


Ve — y ZAv + Vr + Vr 

^ s f ^ a i w p i ^ p 


1 


32 $ 2 foi • va)(v9 ■ vi) 

" aV.Fl F4)\F^ FA) 


V.s ^ "a + "a + "a 


1/2 


64 i^(pi -p 4 )(p 2 -p 3 ) 


v s — > v e + e + + e~ 


1 


128 ^[y|(pi-p 3 )(p 2 -P4)+ 






0h(Pi • P4)(P2 • Ps) + 9l9r™\(V\ ■ pa)] 


^ ->■ ^(r) + e + + e~ 


1 


128 ^ (r) [^I(pi-P3)(p 2 -P4)+ 






0b(Pi • P4)(P2 • Ps) + 9l9r™\{V\ • p 2 )] 



Table 4: Squared matrix elements for decays of sterile neutrinos iV. Here 5 1 is the 
symmetrization factor; a, (3 = e, fi,r; a ^ (3. $ a is the squared mixing angle of 
sterile neutrino with z/ Q . The results are both for Majorana and Dirac neutrinos, for 
details see text. 



^ l^a) — J^VaiWi) (see e.g. [50J for reviews): 
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cos 0y and sy = sin #y are functions of the active-active neutrino mixing 



here 
angles 0^. 

Exact treatment of active neutrino oscillation in the early Universe is a difficult 
task (see e.g. [53l|5TJ[52]) Characteristic timescale of oscillation between i and j mass 
eigen-states for a neutrino with energy E is [50] 



Ti, 



I 2 

\mf 



8.3 x 10" b s 



£ 10" 3 eV 2 



MeV mf 



(20) 



Average energy of relativistic Fermi particles in equilibrium is (E) = 3.15T [33J. 
Applying this relation to active neutrinos and using their measured mass differences 



m 2 . — m\ « 7.6 x 10 5 eV 



Tl2 



1.0 x 10" 3 sec 



|m| — m 
T 

3 MeV ' 



2.5 x 10~ 3 eV , we obtain 



3.1 x 10" 5 sec 



T 



3 MeV 



(21) 



provided that influence of the surrounding environment on neutrino propagation 
is neglected. One sees therefore that about the moment active neutrino decouples 
T ~ 3 MeV typical oscillation timescales are much smaller than the Hubble expansion 
time given by Eq. (J2]) 



th 



15 



Wg m G N T4 
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(22) 
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Here g* ~ 11 (at T ~ MeV) [33] is the so-called number of relativistic species that 
enters energy-temperature relation p = - |* T . Therefore, active neutrinos oscillate 
many times between the subsequent reactions involving them. In quantitative terms 
it means that probabilities P a p to transform from flavour a to flavour f3 are oscillating 
functions of time. In realistic situation neutrinos do not have a definite momentum 
but are created in wave packets that are superpositions of states which have one. 
Since oscillation periods are momentum-dependent according to Eq. f )20|) . each state 
in the superposition will have his own period. Therefore after sufficiently many 
periods initial phases characterizing superposition will change, and there is no reason 
for the phase changes to be correlated with each other. So the decoherence of states is 
what happens. This phenomenon can be described effectively by averaging transition 
probabilities P a p over time. Resulting expressions are [50] 

Pee = 1 — 77 (sin 2# i3 + cos 4 #13 sin 2 20i 2 ) (23a) 

P e)1 = P„e = l cos 2 # 13 sin 2 2# i2 (23b) 

P eT = P Te = sin 2 # 13 cos 2 # 13 sin 2 29 12 ^j (23c) 

P w = l-i S in 2 2# 12 (23d) 

P^r = P T» = \ sin 2 #13 sin 2 29 l2 (23e) 

P TT = 1- sin 2 #13 ^2 cos 2 #13 + ^ sin 2 # 13 sin 2 2# 12 ^ (23f ) 

To understand what happens with a neutrino, consider example of electron- 
neutrino created in electron-positron annihilation. At the production time this par- 
ticle has probability 1 to oscillate into u e and zero for other final state. After long 
enough time for many oscillations to happen and before the time when a collision 
with other particle becomes quite probable, the decoherence comes into play. So now 
we may find the v e with probability P ee , with probability P etl and v T with P er . 
The production rate of the initial specimen per unit time is proportional to collision 
integral I e , according to the Boltzmann equation (|7|). But the actual number of 
produced electron neutrinos is actually reduced by factor P ee . And even if (consider 
this hypothetical situation) muon neutrino does not interact with plasma, it will be 
anyway produced, at rate P ef il e - Generalization to other neutrino flavours leads us 
to conclusion that the modified Boltzmann equation 

describes neutrino dynamics correctly (that is not the case for the initial equation 
©)• 
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However, in dense medium oscillations proceed differently due to considerable 
effects of the plasma on properties of a single particle. Still, the phenomenon can 
be described by the formalism of the PMNS matrix. The difference is that mixing 
parameters together with masses now depend on properties of the environment. In 
case of plasma close to equilibrium with no non-trivial conserving charges present 
the parameter describing it is the temperature. So the parameters of the PMNS 
become temperature- dependent. In language of the effective Hamiltonian approach 
the system of three neutrinos is described by the addition of medium potential AHm 
to the Hamiltonian Hy of the system in vacuum [55J 

H M = H V + AH M , H Y = ^V*dmg(mlmlmj)V^ (25) 

where E is the neutrino energy. Diagonalization of the total propagation Hamiltonian 
Hm gives effective masses and mixings. 

The medium potential comprises effects of neutrino interactions. Since neutri- 
nos take part only in charged- and neutral-current interactions, matter potential has 
two terms AHcc and AH^c, respectively. All neutrinos couple to neutral currents 
identically, so AH^c is proportional to unit matrix. Therefore this term just renor- 
malizes energy, and does not affect oscillations. In contrast, the charged- current term 
is non-diagonal and is present only for v e . The reason is that due to abundance of 
electrons in plasma, v e couples effectively to charged currents, while at temperatures 
below the muon's mass there is no significant contribution of muons and tau-leptons 
to realize coupling of other neutrinos to W boson. 

Explicitly matter potential is [20] 



AHcC = ~^M^ E T ' diag(1 ' °' 0) (26) 



in the flavour neutrino basis (i/ e , u^, v T ). My/ is the mass of the W-boson. 

So far we have dropped sterile neutrinos from consideration. But their mixing 
properties are also altered in hot plasma. Using the approach of the effective Hamil- 
tonian for them, one finds that their effective mixing angles in medium 9 m differ 
from that in vacuum 9y as 



"ji^L „ ^flL „ l0 -n x ( T y ZlOMeVy 

0v M^MI \100lleVj \ M s ) K ' 

for small mixing angles By. Therefore the mixing angle is not altered significantly 
for sterile neutrinos and matter effects are negligible for their dynamics. 
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